Noiseless Quantum Circuits for the Peres Separability Criterion 
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In this Letter we give a method for constructing sets of simple circuits that can determine the 
spectrum of a partially transposed density matrix, without requiring either a tomographically com- 
plete POVM, or the addition of noise to make the spectrum non-negative. These circuits depend 
only on the dimension of the Hilbert space and are otherwise independent of the state. 
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Introduction There has recently been much interest in 
measuring functions of the density matrix directly, with- 
out first performing full state tomography. The non-local 
properties of density matrices have attracted particular 
attention and there has been a series of papers explor- 
ing this possibility [111111111,0]. These typically rely 
on a combination of the Structural Physical Approxima- 
tion (SPA) 7] followed by measuring the spectrum of the 
resulting density operator using the method in (but 
see also II S i). The SPA is a method for modify- 
ing a map which is not a completely positive map (CP 
map) so that it becomes one. This is done by forming 
a convex mixture of the original map with another map 
that projects the state of the system onto the maximally 
mixed state, M{p) = 1/d, for all p. For maximum sensi- 
tivity, we should use as little M{p) as possible. For the 
partial transpose operation on a ci x d dimensional system 
lyi ig) As, this optimal CP-map is given by (1| 

1 

0= ^ Ma<E)Mb+ „ , ^ Ia^Ab. (1) 



d^ + l 



For a more detailed discussion of the properties of this 
map and its implementation as a positive operator- valued 
measure, see |9|. 

This approach was partly inspired by tjie interferom- 
eter circuits in and the analysis in ^0j. The proof 
in ^3 showed that circuits of the form in figure ^ where 




FIG. 1: General form of an interferometer circuit. 

the density matrices undergo the evolution p i~+ [/ pU^ 
conditioned on the state of the control qubit, can mea- 
sure in the modulation of the interference pattern of that 
qubit, 



where v is the visibility, and a is a phase shift. This is sig- 
nificant, because it enables us to characterise the effect of 
maps defined by a unitary map on the kets only, without 
being constrained to perform W on the bra vectors. 

If we make p — p®p®...®p from multiple copies of 
p, an m X m density matrix of interest, we can use this 
circuit to measure the moments of the density matrix, 
Tr(p^) . . . Tr(p™) by choosing U to be the cyclic shift |1| 

y«|^i)|^2)...|^;.)-|^fe)l¥'i).-.|¥'fc-i), (3) 

acting on the kets of the basis only, thus measuring 

m 

rY^p^ky(k) ^ rj^^pk ^ ^ , (4) 

1 = 1 

from which we can determine the eigenspectrum of the 
density matrix P, H ^ . 

Let us begin by considering polynomial local-unitary- 
invariants of pure states. These can be constructed as 
follows. Begin by writing 
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where the tensor indices denote the various subsystems 
in the usual way. Now contract every "downstairs" (or 
"input" ) index with one of the "upstairs" (or "output" ) 
indices labelling the same subsystem, (but not necessar- 
ily on a neighbouring term) until no free indices remain. 
These functions are therefore invariant under local basis 
changes by construction. 

In 1121] Rains proved that there was a one-to-one cor- 
respondence between the set of all local polynomial in- 
variants and sets of permutations that permute the down- 
stairs indices corresponding to each subsystem with other 
downstairs indices corresponding to the same subsystem. 
(Or equivalently, perform the inverse permutation on the 
upstairs indices.) Rains' work was motivated by the need 
to find new shadow enumerators of quantum codes and it 
was not until |l3| that these invariants were given a phys- 
ical interpretation in terms of the expectation values of 
operators. In 0| the methods of [J were extended to 
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include arbitrary permutations, thus recovering Rains' 
construction for the full set of polynomial invariants. 

The preceding papers seem to share the assumption 
that the only way to measure the effect of a map on p 
is to cause the density matrix to actually undergo that 
map as a physical evolution: the map must be physically 
implementable. For maps that are not in this class one 
must therefore use the SPA j7|. (However, we will argue 
later that a more subtle measurement model is implied 
by certain results in some of these papers.) 

This assumption makes intuitive sense and one might 
think that it is therefore "obviously" true; except, like 
some other "obviously true" assertions about quantum 
mechanics, it is not. To illustrate this, we will give a 
method for constructing a set of unitary circuits that can 
measure the spectrum of p^^ for any p, without adding 
noise of any kind. The structure of these circuits de- 
pends only on the dimension of p and does not require a 
tomographically complete POVM at a ny s tage 0. 

Measuring the Kempe Invariant In |l4| a method for 
constructing circuits to measure polynomial invariants 
was presented. Although the authors do not mention 
this fact, this class of circuits includes a set that can 
measure the effects of some non-positive maps obtained 
by extending maps that are positive, but not CP maps 
to act on larger spaces. 

Consider polynomial invariants of pure states of three 
qubits. A 6*'' order invariant (in ip) that is algebraically 
independent of all the lower-order invariants was discov- 
ered by Kempe and can be written (with a little 
rearranging) as 

h = Vkipip VimqW VjnrV ■ (6) 

This index-based notation is rather opaque. Figure 
is the same invariant drawn as a summation diagram: 
the linked nodes represent summations; each upstairs- 
downstairs pair of nodes after each copy of p denotes 
the indices for one local subsystem (qubit), just as with 
tensor index notation. The triangles represent indices la- 
belling qubit 1; the squares, those for qubit 2 and the 
open circles correspond to qubit 3. 




FIG. 2: The summation diagram for the Kempe invariant. 

We can always choose to order the copies of p in the in- 
variant so that one subsystem is traced out; this is just 
another manifestation of a theorem due to Schrodinger 
161 llTl Il8| , which is also known as the GH JW theorem 
19t l2(i l2lj . The traced out party is represented by the 
open circles. This way of writing invariants highlights 
the Rains correspondence between invariants and permu- 
tations of the indices for each subsystem Consider 




FIG. 3: Summation diagram for another 3-qubit invariant, 
the following three qubit invariant: 

This is a three qubit invariant but we are interested in 
bipartite entanglement here. If we ignore the third qubit 
altogether, we can see that this invariant is Tr((pi2)'^), 
where the subscripts 1, 2 are no longer indices, but now 
refer to qubits 1 and 2 and the two longest lines represent 
the final matrix trace. We are free to drop the third qubit 
because the indices for particle 3 are always traced out in 
equation ^ and their links in the summation diagram 
(Fig. I^J are never exchanged between copies of p. 

We can now connect figure [3 to the partial transpose 
map. One way to obtain|21from|31is to take each upstairs- 
downstairs pair of indices for particle 2 and swap the up- 
stairs index with the downstairs one. This performs the 
partial transpose operation on each copy of p. Therefore 

/6=Tr((pS)'), (7) 

where the symbol denotes the partial transpose with 
respect to qubit 2. In other words, equation © is the 
third moment of the partial transpose of p. 

We can measure the Kempe invariant because the sum- 
mation in Fig.|5]can be obtained from that in Fig.|21in two 
equivalent ways. The other way is via the Rains construc- 
tion, using a cyclic permutation acting on the downstairs 
indices for qubit 2 to reproduce the "contraflow" summa- 
tion pattern in Fig. |21 As all the indices are summed out 
in upstairs-downstairs pairs, their labels have no physi- 
cal significance and so these two operations produce the 
same invariant. Any invariant of the form in Q with no 
free indices and all summations on upstairs-downstairs 
pairs can be generated using Rains' construction [l^ . 

A circuit can be obtained for this invariant by using the 
construction in [l^ . This is done by rotating the sum- 
mation diagram in figureElby 90° clockwise and omitting 
the p symbols (see Fig. 0]) The downstairs indices now 
correspond to input rails and the upstairs indices to out- 
put rails. 

A > 

O ' \ O 

*/ 

o o 

FIG. 4: The Kempe invariant, in the notation from [l^ . 

To construct the circuit, insert the wiring in Fig. ^ as the 
controlled-?/ in Fig. 
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Measuring the spectrum of p^^ It can be seen that we 
can also write the invariants Tr((/9j''|)^) and Tr((p]'2^)'*) 
in a similar form, so we can construct circuits to evaluate 
these moments. We can therefore determine the eigen- 
spectrum of using the recipe given in 0,0] (or alter- 
natively, reconstructing the characteristic polynomial for 
from the values of Tr((p^2^)'^)). Thus we do not need 
to use the SPA to determine the spectrum of p^'^ . 
The circuits for the second, third and fourth moments 
are shown in Figs. I5I7I respectively, with the permuta- 
tions decomposed into controUed-SWAP (CSWAP) gates. 
These are the only circuits required for two qubit states; 
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FIG. 5: The circuit for TvUpJ^f'' 
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FIG. 6: The circuit for Tr((p^2')^) 




FIG. 7; The circuit for Tr((p^2' )'')• 

we do not need a circuit for Tr(p^| ) because the partial 
transpose is a trace preserving map, so the first moment 
is already known to be 1. In fact, the only unitary circuit 
we can write down that is first order in p, or indeed 
is the circuit for the trace norm. This is also the only 
polynomial invariant that is quadratic in tp. 



Discussion We have exhibited a class of circuits that 
can measure the spectrum of p"^^ , without needing to per- 
form a tomographically complete POVM or add noise. 
While we have only given the circuits for the two qubit 
case, the construction generalizes to bipartite systems of 
any finite dimension in the obvious way; all that is re- 
quired is two cyclic permutations of opposite handedness 
acting on the rails for each subsystem. Fiurasek |^ has 
shown that the SPA for the partial transpose for two 
qubits is a tomographically complete CP map. It is an 
open problem whether this is generally true for the par- 
tial transpose on larger systems, but the method given 
in this paper never needs a tomographically complete 
POVM to characterize the effects of the partial trans- 
pose on p, regardless of its dimension. These circuits are 
surprisingly simple: if we can implement an eigenspec- 
trum measurement by the method in P, Q then we can 
measure the spectrum of p"^^ with a set of circuits that 
uses the same set of primitives. The complexity of the 
networks also scales the same way as those for eigenspec- 
trum circuits; for a general di x d2 dimensional mixed 
bipartite state, the circuits consume the same number of 
copies per run, ^d^d"^ + \d1d2 ~ 1. 

What may be more interesting about these circuits is 
that they disprove our intuition that the only characteri- 
zable map is a physically implementahle map. This state- 
ment is still true if we only have access to one copy of the 
state at a time, as this method for measuring the effects 
of maps needs multiple copies in order to work; since we 
cannot implement the map directly, we must characterise 
it via its moments. It is also worth noting that having 
access to only two copies at a time is not much use either; 
the circuit for Tr((y0^2^)^) is identical to that for Tr(p^2)- 
Thus the partial transpose may also preserve Tr(yO^). 

If we can manipulate at least did2 copies at once, we 
can perform the Peres test This is a necessary and 
sufficient condition for separability for two qubit states, 
although it is not sufficient to test for separability in gen- 
eral [231 • The entire minimal generating set for the ring 
of polynomial invariants described by Rains 0, Q| is 
necessary to completely separate the orbits of the state 
under local unitary operations in general. It is also suf- 
ficient 0, . However, there is no known upper 
bound on the order of the invariants needed to construct 
a complete minimal generating set. One way to see this 
is to notice that this paper has only examined invariants 
whose summation pattern has period 1 in the number 
of copies of the density matrix it passes through before 
it repeats. There is no a priori reason why there might 
not be invariants corresponding to longer-period summa- 
tion patterns that cannot be written as products of terms 
with period 1 and thus could have a much larger order. 

This is not the only reason why this problem is hard, 
but it is enough to show that a general density opera- 
tor may need circuits in more than did2 copies in order 
to characterize its non-local properties completely. We 
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cannot say how many copies we would need to be able 
to process at once without first finding the minimal gen- 
erating set of invariants for the system. This is a hard 
problem and all we can say in general is that this number 
is finite [13 . 

If we assume that we can use unitary circuits that can 
manipulate any finite number of copies at once, the class 
of measurable maps is considerably larger than the set 
of all unitary maps. Consider a density matrix that we 
have the option of pre-processing according to 

for any unitary V. If our interferometer circuit contains 
a controlled- [/ of our choice, we can now measure any 

p®" ^ T/Fp®"F'f (9) 

where W = UV, which is a much larger class of oper- 
ations than is accessible by equation JHJ alone. What 
is not immediately obvious is that this class of oper- 
ations is qualitatively different from the usual unitary 
operations. This is because we now have the ability to 
operate on the ket basis of the density matrix without 
having to perform the adjoint operation on the bra basis. 
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